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Extracting Galaxy Cluster Gas Inhomogeneity from X-ray Surface 
Brightness: A Statistical Approach and Application to Abell 3667 
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ABSTRACT 

Our previous analysis indicates that small-scale fluctuations in the intracluster 
medium (ICM) from cosmological hydrodynamic simulations follow the lognormal 
probability density function. In order to test the lognormal nature of the ICM 



^ ! directly against X-ray observations of galaxy clusters, we develop a method of 

- ' extracting statistical information about the three-dimensional properties of the 



fluctuations from the two-dimensional X-ray surface brightness. 



O 

c/3 , We flrst create a set of synthetic clusters with lognormal fluctuations around 

their mean proflle given by spherical isothermal (3 models, later considering poly- 
tropic temperature proflles as well. Performing mock observations of these syn- 
^ ! thetic clusters, we flnd that the resulting X-ray surface brightness fluctuations 

^ i also follow the lognormal distribution fairly well. Systematic analysis of the 

^ [ synthetic clusters provides an empirical relation between the three-dimensional 

►v^ I density fluctuations and the two-dimensional X-ray surface brightness. 

O I We analyze Chandra observations of the galaxy cluster Abell 3667, and flnd 

OO 

Q I that its X-ray surface brightness fluctuations follow the lognormal distribution. 

While the lognormal model was originally motivated by cosmological hydrody- 
'k>( I namic simulations, this is the flrst observational conflrmation of the lognormal 

H I signature in a real cluster. 

Finally we check the synthetic cluster results against clusters from cosmologi- 
cal hydrodynamic simulations. As a result of the complex structure exhibited by 
simulated clusters, the empirical relation between the two- and three-dimensional 
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fluctuation properties calibrated with synthetic clusters when applied to simu- 
lated clusters shows large scatter. Nevertheless we are able to reproduce the true 
value of the fluctuation amplitude of simulated clusters within a factor of two 
from their two-dimensional X-ray surface brightness alone. 

Our current methodology combined with existing observational data is useful 
in describing and inferring the statistical properties of the three dimensional 
inhomogeneity in galaxy clusters. 



Subject headings: galaxies: clusters: general 
mology: observations 



X-rays: galaxies: clusters - cos- 



Introduction 



Galaxy clusters have been one of the most important probes of cosmology (e.g., Bartlett 
& Silk 1994; Eke, Cole, & Frenk 1996; Viana & Liddle 1996; Kitayama & Suto 1996, 1997; 
Kitayama, Sasaki, & Suto 1998; Holder et al. 2000; Haiman, Mohr, & Holder 2001; 
Majumdar & Mohr 2004). In the context of dark energy surveys, which attracts much of the 
attention of the cosmology and particle physics communities, galaxy cluster surveys are also 
unique in that they most direct ly probe the growth of structure rather than relying solely 



on distance measurements (e.g.. lAlbrecht et al. 1120061 ). In order to capitalize on this, galaxy 



cluster surv eys, in particular those utilizing the Sunyaev-ZeFdovich effect (for reviews see 



for ex ample, ICarlstrom et al.ll2002l : iBirkinshaw Ill999l : iRephaeli Ill995l : ISunyaev fc Zel'dovich 



1980l ). are currently operating and many more are planned in the near future. However, 
one of the biggest challenges in interpreting these surveys is relating physical quantities 
of galaxy clusters, namely mass, to observable ones. In particular, these mass-observable 
relations may be sensitive to the inherent complex structure of clusters. Therefore we must 
better understand galaxy clusters to utilize fully the potential of galaxy cluster surveys in 
constraining cosmological parameters. 

Recent observations of galaxy clusters have revealed a rich variety of structural complex- 
ity. Recent X-ray satellites with their improved angular resolution, collecting area, and si- 
mult aneou5_spectra]_measuiwie^ have unveiled complex ter nperature structur e 
fe.g.. iMarkevitch et al. II2OOOI: iFurusho et al. Il200ll ). shock front s fe.g.. I Jones et al.l 20021 ). 



cold fronts (e.g., IMarkevitch et al. Il2000l ). and X-ray holes (e.g., iFabian et al. II2OO2I ). Im 



proved observational strategies and analysis methods of lensing observations of galaxy clus- 
ters s how that the rnass d istribution, as opposed to just the gas, is often complicated as well 
(e.g.. lBradac et al. II2006I ). Both X-ray and lensin g observations of g alaxy clusters reveal that 



clusters are frequently undergoing mergers (e.g., iBriel et al.l 120041 ). With such various and 
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sundry structural complexities may galaxy clusters reliably be used as cosmological probes 
? 

The complex structure seen in galaxy clusters motivates our investigation of the intra- 
cluster medium (ICM) inhomogeneity. We note, however, that we take a statistical approach 
to modeling the inhomogeneity rather than directly modeling such complex phenomena as 
shocks, cold fronts, etc. Motivated by results from cosmological hydrodynamic simulations 
we explore the ramifications of a lognormal model of the inhomogeneity of the ICM. This 
model was first proposed (Kawahara et al. 2007, hereafter Paper I; Kawahara et al. 2008), in 
this context, to explain the discrepancies between emission weighted and spectro s copic tem- 



perat ure estimates from galaxy clusters (JMazzotta et al.l l2004l : iRasia et al.l l2005l : IVikhlinin 



20061 ) . They found that local inhomogeneities of the ICM play an essential role in producing 



the systematic bias between spectroscopic and emission weighted temperatures. 

Thus far, the lognormal model has been motivated by and applied only to clusters from 
cosmological hydrodynamic simulations. Therefore it is crucial to see if inhomogeneities in 
real galaxy clusters also show the lognormal signature. In reality, this is not a straightfor- 
ward task since one can observe clusters in X-rays only through their projection over the 
line of sight. Thus we develop a method of extracting statistical information of the three- 
dimensional properties of fluctuations from the two-dimensional X-ray surface brightness. 

The rest of the paper is organized as follows. We first summarize the log-normal model 
in ^ We create synthetic clusters to explore the relationship between the intrinsic cluster 
inhomogeneity and X-ray observables in ^ In §1] we apply our methodology to Chandra 
observations of the galaxy cluster Abell 3667, and then attempt to quantify the nature 
of cluster inhomogeneity. We also compare our synthetic cluster results with cosmological 
hydrodynamic simulations in §3 Finally, we summarize our results in §0 Throughout the 
paper the Hubble constant is parameterized by h in the usual way, Hq = 100 h km s~^ 
Mpc-^ 



2. Model of the ICM Inhomogeneity 

2.1. Lognormal Distribution 

In order to characterize the inhomogeneity of the ICM, we define the density and tem- 
perature fluctuations as the ratios 6n = n{r)/n{r) and 5t = T(r)/T(r), where n(r) and T(r) 
are the local density and temperature at radius r, and n{r) and T{r) are the angular average 
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profiles defined by 

If 

nir) = — / n (r, 6, (p) sin 9 d9 d(j) (1) 

An J 

T{r) = -^ [t (r, 9, (J)) sine dOdcj) (2) 

where 6 and (f) are polar and azimuthal angles, respectively. Analysis of hydrodynamical 
simulations (Paper I) found that 5„ and 6t are approximately independent and follow the 
radially independent lognormal probability density function (PDF) given by 

,21 



Pi^x] o-LN,x) dSx = —1= exp 






0-r. 



-'X 



where x denotes n or T, Sx = x{r)/x{r), and ctln.x is the standard deviation of the logarithm 
of density or temperature. 

To construct the two-dimensional surface brightness profile from the three-dimensional 
density and temperature distribution, we also need the properties of the power spectra of 
the density and temperature fluctuations. We adopt statistically isotropic fluctuations with 
a power-law type power spectrum for both the density fluctuations Pn{k) oc fc"" and the 
temperature fluctuations Pxik) oc k""^. These assumptions are based on the results of the 
cosmological hydrodynamic simulations described in §2.21 



We use this model to generate synthetic clusters to explore the relationship between 
the three-dimensional inhomogeneity in the ICM and the two-dimensional X-ray surface 
brightness. 



2.2. Cosmological Hydrodynamic Simulated Clusters 

When one considers the projection of galaxy clusters to two dimensions for mock X-ray 
observations, the power spectrum of the fluctuations is important in addition to the PDF 
of the inhomogeneity. Here, we once again turn to simulations to investigate the power 
spectrum of the fluctuations. 

We extract the six mas sive clusters from co smological hydrodynamic simulations of the 



local universe performed by iDolag et al.l (120051 ). The simulations utilize the smoothed par- 



ticle hydrodynamic (SPH) method, and assume a flat A CDM universe with flrn = 0.3, flh = 
0.04, as = 0.9, and a dimensionless Hubble parameter h = 0.7. The number of dark matter 
and SPH particles is ~ 20 million each within a high-resolution sphere of radius ~ 110 Mpc, 
which is embedded in a periodic box ~ 343 Mpc on a side that is filled with nearly seven 
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million low-resolution dark matter particles. The simulation is designed to reproduce the 
matter distribution of the local universe by adopting the initial conditions based on the IRAS 
galaxy distribution smoothed over a scale of 4.9/i~^Mpc. Thus, the six massive clusters are 
identified as Coma, Perseus, Virgo, Centaurus, A3627, and Hydra. A cubic region with 6 
h~^ Mpc on a side centered on each simulated cluster is extracted and divided into 512^ cells. 
The density and temperature of each mesh point are calculated from SPH particles using the 
B-spline smoothing kernel. A detailed description of this procedure is given in Paper I. The 
distance between two adjacent grid points is given by t/grid = 6/i^^Mpc/512 ~ 12h^^ kpc, 
which is comparable to the gravitational force resolution (14 kpc) and the inter-particle sep- 
aration reached by SPH particles in the dense centers of clusters. Therefore, the (maximum) 
resolution is dgj-i^/rc ~ 0.1 assuming Tc ~ 100 kpc. This is about one order of magnitude 
worse than that of both the synthetic clusters (§ [3]) and the observational data (§ Hj). 

For each simulated cluster, we compute the radially averaged density and temperature 
profiles, n{r) and T{r), respectively (Eqn. [1] and [2]), and use them to compute the density 
and temperature fluctuations Sn = n/n and 6t = T/T at each grid point. We extract 128'^ 
cells of 6n and St around the center of a simulated cluster and compute the power spectrum. 
The distance from the center to the corner of the 128^ cells is ~ 1.3 h~^ Mpc which is 
approximately equal to the virial radius of the simulated clusters (r2oo = 1.0-1.6 h~^ Mpc). 
The virial radius, r2oo, is the radius within which the mean interior density is 200 times that 
of the critical density. 

Figure [H shows the power spectra for each simulated cluster for both 6n (upper panel) 
and St (lower panel). In each panel a simple power law, P{k) oc k~^ (dotted line), is 
also plotted for comparison. The power spectra for both the density and temperature are 
relatively well approximated by a single power law. We therefore adopt a power-law spectral 
model for the density and temperature fluctuations for the synthetic cluster analysis. 



3. Synthetic Clusters 

Cosmological hydrodynamic simulations provide a useful test-bed for exploring cluster 
structure. Simulated clusters exhibit complex density and temperature structure akin to 
that of real galaxy clusters. The resolution of our current simulations, however, is limited, 
especially when compared to the resolution available from current generation X-ray satellites. 
In addition, we need to systematically survey the parameter space of (TLN,n and a„ in order 
to relate the X-ray surface brightness fluctuations to the density fluctuations. Thus we 
create a set of synthetic clusters at higher resolution that have lognormal fluctuations around 
their mean profile. Analysis of mock observations of these synthetic clusters enables us to 
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Fig. 1. — The power spectra of 5„ (upper) and 5t (lower) of the six simulated clusters. 
Dashed lines indicate P{k) oc k~^. 

investigate the relation between the X-ray surface brightness and the statistical properties 
of the three-dimensional density fluctuations, namely o"LN,n and a„. 



3.1. Method 



3.1.1. Synthetic Cluster Generation 

The three-dimensional synthetic clusters will be projected to two dimensions when con- 
sidering the X-ray surface brightness. In order to incorporate a power-law type power spec- 
trum of spatial fluctuations into the synthetic clust ers, we follow a similar methodology 



as th at of several studies of the interstellar medium (jElmegreen I l2002l : iFischera &: Dopita 



2004 ). First a Gaussian random field with a power-law power spectrum is constructed and 
that field is mapped into a lognormal field. Therefore, our assumption for the power spec- 
trum is adopted for the Gaussian field q as opposed to 5n. However, we will verify that the 
ensemble average of the power spectra of q and 5n {Pq{k) oc /c"' and P„(/i;) oc A;"") have 
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almost the same power-law indices, a^ ~ q;„. 

We generate the lognormal density fluctuation field as follows. We first generate the 
real random fields, a(k) and 6(k), in fc-space, whose distribution functions obey 



p{a)da 



1 



V^^m 



exp 



fik) 



da, p{b)db 



1 



V^^m 



exp 



m 



db. 



(4) 



where f{k) = Ak°''' . Then we compute q{r), the Fourier transform of a complex field 
g(k) = a(k) +ib{k). With the additional conditions a(k) = a(— k) and 6(k) = — 6(— k), g(r) 
becomes a real Gaussian random field, and its power spectrum, Pq{k), is equal to the input 
function f{k) = Ak"""^. The amplitude A is related to the variance of the Gaussian random 
field: 



^s = 4vr 



Pf{k)dk, 



(5) 



where k^[^ and fcmax denote the minimum and maximum value of the wavenumber. Finally 
the lognormal deviate, Sxir), is obtained from the Gaussian deviate, q{r), using the relation 



exp 



0"LN,a: 



a 



(Tr, 



-q{r) - 



LN,x 



(6) 



where ctln, x is the standard deviation of the lognormal field. 



We construct synthetic clusters with average density given by the (3 model and (5„ drawn 
from a lognormal distribution taking into acc ount the power-law type power spectru m of 
spatial fiuctuations. The /? model is given by ( jCavaliere fc Fusco-Femiano Ill976l . 119781 ) 

-3/3/2 



n(r) = no 



1 + 



(7) 



where Uq is the central electron number density, Vc is the core radius, and /3 specifies a 
power-law index. For simplicity, we first adopt a fiducial value of /? = 2/3, and assume 
isothermality for the synthetic clusters. Later, we examine the effects of varying /3 (§ 13.3.11) 
and of temperature structure using a polytropic temperature profile (§ 13.3.21) . 

The density at an arbitrary point is given by 

n(r) = (5„n(r). (8) 

The X-ray surface brightness profile is obtained by projecting the three-dimensional syn- 
thetic cluster down to two dimensions. For the isothermal case the projected X-ray surface 
brightness profile is 

5x(R) oc f[n{v)fdl, (9) 
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< 
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Fig. 2. — The change of the power-law spectral index of the density («„) and density squared 
fields (a„„) compared to that of the Gaussian field (og). Solid and dashed lines indicate 
«n/ttq — 1 (density) and ann/c(q — 1 (density squared), respectively. Each symbol indicates 
a different value of o"LN,n (cross, square, and triangle correspond to ctln,™ = 0.1, 0.3, and 
0.5, respectively.) The power-law index of the density field is very close (< 3%) to that of 
the Gaussian field used to generate the lognormal distribution and that of the square of the 
density is within ~ 13% for larger values of ULN.n and < 5% for smaller values (crLN,n ^ 0.3). 



where R indicates the position on the projected plane and / is the projection of r onto the 
line of sight direction. 

Performing the procedure described above, we set up a cubic mesh of n{r) in which our 
three-dimensional synthetic cluster is located with A^grid = 512 grid points along each axis. 
We choose the box size Lbox = 10 r^., which results in the distance between two adjacent grid 
points being dgrid = lOrc/A^grid ~ 0.02 Tc. 

We fit the power spectrum of the 5„ field by a power-law spectrum so that P„(k) oc A;"". 
We also fit the power spectrum of the square density field, 6nn = v? j {v?) = 5„^ exp (— o"ln „) 
(Appendix B), by the power-law Pnn{k) oc k""", relevant to X-ray surface brightness since 
Sx oc f di in?. Throughout this paper, the notation (x) is used to denote the ensemble 
average of quantity x over many clusters. 

Figure [2] shows the change of the power-law spectral index of the density (q;„) and 
density squared fields (a„ri) compared to that of the Gaussian field (a,). The change in the 
power-law index for the density and density squared distributions compared to the initial 
Gaussian field are small (<3 % and < 13%, respectivel y), and therefore, a 



~ an ~ a-r, 



consistent with the results of iFischera &: Dopita I (120041 ) 
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3.1.2. X-ray Surface Brightness 

To quantify the relationship between the inhomogeneity of the density and the X-ray 
surface brightness, Sx, we introduce the X-ray surface brightness fluctuation from the average 
radial surface brightness profile S'x(-R) 



5sx(R) 



gx(R) 
SxiR)' 



(10) 



where R = |R|. We define the average profile Sx{R) for an individual cluster by fitting the 
projected synthetic clusters to an isothermal /3 model 



Sx{R) — Sx,{ 



R 

rc.x 



-1 -3/3X+1/2 



;ii) 



where Sx,o is the central X-ray surface brightness, r^x is the core radius, and /?x specifies 
the power-law index for the X-ray surface brightness distribution. These three parameters 
are derived from a model fit to each synthetic cluster. It is important to emphasize that the 
average in equation (ITT!) is defined for an individual cluster. We note that if we adopt directly 
the average X-ray surface brightness profile instead of a /3 model fit (Eqn. [H]), the results 
are unchanged. This is because the radial profile is well approximated by the /3 model 
for the synthetic clusters. However, for observations of real galaxy clusters, the j3 model 
approximation might break down and one should instead use an average of 5'x(R) directly 
in such cases. In §3.2, we will investigate the relation between the standard deviation of the 
X-ray surface brightness fluctuations, uln, Sx, and that of the intrinsic density fluctuations, 

0"LN,n- 

Here, we consider the relation of ctln, Sx and o"LN,n for the ensemble average of clusters 
assuming they all obey the /3 model with the same /3, re, aq and crLN,n: 



{Sx){R) = (Sx(|R|)) 
(5x(R)) oc e"£N,n / n^dl 



(12) 
(13) 



where the exponential term of the right hand side of equation (IT3l) comes from the second 
moment of the lognormal distribution (Paper I). Although the ensemble average is not an 
observable quantity, we can describe an analytical prediction of crLN,Sx(-R) assuming the 
isothermal P model (Appendix |A]) . In addition, one expects that 5'x ~ (5'x) if there is a 
large enough volume compared with the size of fluctuations when calculating 5'x. In other 
words, the spatial average approaches the ensemble average. For these reasons, it is useful 
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to consider the ensemble average. Using equations flT2l) and flT3l) . we define the ensemble 
average of fluctuations in the X-ray surface brightness as 

We note that the distribution of the square of density fluctuations, which is proportional 
to the local emissivity in the isothermal case, is also distributed according to the lognormal 
function with a lognormal standard deviation of 2crLN, n if the density fluctuations follow the 
lognormal distribution with standard deviation ctln.h (Appendix [B|) . 



3.2. Statistical Analysis of the Synthetic Clusters 

Here, we investigate the distribution of (5sx of the synthetic clusters and relate quantities 
obtainable from observations, ctln, Sx and asx, to that of the underlying density, o"LN,n and 



3.2.1. Lognormal nature and the relation between ctln, Sx CLnd crLN,n 

We investigate the distribution of (5sx as a function of radial distance R from the cluster 
center. We first divide the ^sx field into shells of thickness 0.5 Tc. The distributions of ^sx 
within each shell, p{Ssx',R), averaged over 256 synthetic clusters are shown in Figure [3] for 
various values of ctg. We find that ^sx also approximately follows the lognormal distribution. 
The standard deviation of the logarithm of 6sx versus radius, crLN,Sx(-R), constructed from 
the averaged shells is displayed in Figure HI Two values of cTLN.n are plotted, 0.1 and 0.5, in 
addition to using the average profile defined by both the /3 model (Eq. [TT]; solid) and that 
for the ensemble (Eq. [13]; dotted). The analytic prediction (Eq. jA13] : dashed) and the 
case including the temperature structure ( §3.3.2t dot-dashed) are also plotted. At large R, 
c"LN,Sx,ens(-R) IS approximately crLN,Sx(-R) because the spatial average tends to the ensemble 
average due to the large volume used for averaging. However, the agreement is poor near the 
center, where the ensemble average is not a good approximation. Although only one value 
for aq is shown, similar results are obtained for other values. 

Figures [3] and H] indicate that the probability density function is weakly dependent on 
the projected radius R. This radial dependence is caused mainly by two competing effects. 
Consider first the case where the typical nonlinear scale of fluctuations is much smaller than 
the size of the cluster itself (shallow spectrum). As equation (Al) indicates, the surface 
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Fig. 3. — The probability distribution of the ensemble-averaged distribution of 6sx illustrating 
the radial dependence. The distributions in shells of thickness 0.5 Tc are shown. Each color 
indicates a different radial interval: R < 1.5 r^. (red), 1.5 r^ < R < 3.5 r^ (black), and 
i?>3.5rc (blue). 



brightness at R is given by 



Sx{R) oc / 5, 



1 + 



f 



r2 + i?2 



-1-3/3 



dl. 



(15) 



This implies that the mean value of Sx{R) is effectively determined by the integration over 
the line of sight weighted towards the cluster center, roughly between — -^/r^ + i?^ g^^^ 
+ \/^c + -^^- This is also true for the variance of Sx{R)- Since the effective number of 
independent cells contributing to the variance of Sx{R) is smaller at smaller projected radii, 
c"LN,Sx slightly increases for smaller R. This explains the behavior of the shallow spectra 
results for aq = —2 and —2.5 in Figure [3l On the contrary, if the typical nonlinear scale 
of fluctuations is comparable to or even larger than the cluster size (steep spectrum), the 
sampling at the central region significantly underestimates the real variance. So the ctln.Sx 
should increase toward the outer region. This is seen in Figure 3 for the steeper spectra. 
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Fig. 4. — The radial dependence of the standard deviations of the logarithm of X-ray surface 
brightness, o"ln, Sx- Two values of crLN,n are plotted, 0.1 and 0.5, as indicated in the figure. 
Solid and dotted lines show o"ln,Sx(-R) calculated using the average profile defined by the 
(3 model (Eq. [11]) and the ensemble average (Eq. [13]), respectively. Dashed fines show 
the analytical prediction (Eq. jA13] ) . Dash-dotted lines indicate the case including the 
temperature structure. Although we show results only for a single power-law index, aq = 
—3.0, similar results are obtained in other cases. 



ag = —3.5 and —4. 

Note the first effect is very small and the second effect becomes significant only when 
aq < —3. The cosmological hydrodynamic simulations imply that the typical value of a^ is 
—3. Therefore we neglect the radial dependence of the 6sx field in the following analysis. 

From actual observations, we obtain the 5sx niap for an individual cluster, not the 
ensemble average. Therefore, we evaluate the distributions of (5sx in individual synthetic 
clusters. Figure [5] shows the PDF for five individual synthetic clusters (solid) along with 
the best-fit lognormal distributions (dashed). We neglect the radial dependence and use the 
distribution for the whole cluster within a diameter of Lbox = 10 Tc. Each color represents a 
different individual synthetic cluster and each panel shows a different value of the power-law 
index of the Gaussian field, a^, with values between -2 and -4. Even if the analysis is done 
for one cluster, the distribution approximately follows the lognormal distribution. 



The noisy behavior for steeper spectra 



[Oia 



-3.5, —4) in Figure 5 is due to the 



presence of fiuctuations on scales larger than that of the cluster, similar to the discussion 
above for Figure 3. In other words, steeper spectra {aq < —3) have relatively more larger 
scale fiuctuations compared to shallower spectra {aq > —3). Cosmological hydrodynamic 
simulations suggest that a^ ~ 3, placing galaxy clusters in the less noisy regime. We do not 
consider the noisy regime further in this paper. 



-13- 




Fig. 5. — The probability distribution of 5sx for five individual synthetic clusters (solid) along 
with the best-fit lognormal distributions (dashed). Each color shows a different individual 
synthetic cluster. Each panel shows a different value of the power law index of the Gaussian 
field, aq, between —2 and —4 as indicated in each panel. 



The standard deviations of the logarithm of 5sx, ctln.Sx, for the different sets of a, 
(symbols) and fXLN.n (colors) are shown in Figure [61 The relation between o"LN,n and o"ln,Sx 
is approximately linear (right panel) although the proportionality coefficient depends on aq. 
Therefore, we can write 



C^LN.Sx — <5(ttg)crLN,n- 

We find that Q{aq) can be approximated well by the following function 



Cl 



C2 + \aq\ 



(16) 



(17) 



We calculate the average of crLN,Sx/o"LN,n for each aq over three different values of crLN,n 
(o'LN.n = 0.1,0.3, and 0.5). By fitting o"LN,Sx/o"LN,n(«g) using equation (fTTj) . we obtain 
Cl = 2.05 X 10-2 and Cs = 1.53 x lO'^. 
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Fig. 6. — The average of (Tln.Sx from the 256 synthetic cluster sample as functions of a, 
(left) and o"LN,n (right). The left panel also shows the standard deviation of (Tln.Sx from the 
256 synthetic clusters and black, red, and blue represent different values of ctln,™, 0.1, 0.3, 
and 0.5, respectively. In both panels, symbols indicate values of aq (cross, square, triangle, 
asterisk, and circle correspond to a^ = —2, —2.5, —3, —3.5, and —4, respectively). Dashed 



lines show the best-fit approximately linear o"ln, Sx-<7ln, 
each pair of crLN,n, «g- 



relation (Eq. [11] and Eq. [T7]) for 



3.2.2. Spectral Considerations 

Because ctln, Sx is strongly dependent on the power-law index a^, the estimate of aq from 
the 5sx map is crucial for interpreting the value of o"ln,Sx- Because aq is an un-observable 
quantity, we investigate the relationship between the power spectra of 5n and (5sx by fitting 
the power spectrum of (5sx under the assumptions of both statistical isotropy and a power 
law so that Ps^{K) oc K""^^, where K indicates the two-dimensional wave vector. 

Figure [7] shows the power-law index of the X-ray surface brightness, asx, as a function 
of its counterpart Gaussian field, aq. Averages and standard deviations over 256 synthetic 
clusters are shown for three values of the standard deviation of the logarithm of density, crLN,n, 
where crosses, squares, and triangles correspond to ctln,™ of 0.1, 0.3, and 0.5, respectively. 



The dotted line corresponds to the relation asx 



aq and the solid line shows asx = «g + 0.2. 
We find that asx ~ ctg + 0.2 and since aq ~ «„, this implies asx ~ «n + 0.2. This can 
be understood as follows. As we have seen in § 13.11 the difference between q;„ and ann 
is relatively small (< 13% and often < 5%). If one assumes 5sx is the projection of 5nn 
(although this is only strictly true if the average of the surface brightness is defined by the 
ensemble average as Eq. [12]), 5^^ can be described as 
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Fig. 7. — Comparison of the X-ray surface brightness (ctsx) and the input Gaussian field (oq) 
power-law indices. Symbols and error bars indicate the average and the standard deviation, 
respectively, of asx for 256 samples for different sets of aq and crLN,n- Symbols correspond 
to different values of 0"LN,n5 with cross, square, and triangle symbols indicating aLN.n = 0.1, 
0.3, and 0.5, respectively, and the relations asx = Oo and asx 



Uq + 0.2 are also shown 



(dotted and solid lines, respectively). We obtain asx for each individual synthetic cluster by 
fitting Psy^ (K) of an individual cluster under the assumption of both statistical isotropy and 
a power-law (oc K""^^). 



where indicates celestial coordinates and W{&,1) is the window function. If we neglect 
the ©-dependence of the window function and set W{@,1) = W{1), then Ps^(K.) can be 
written as 

Ps^{K) = ^Jdk Pnn(k) \W{k)\^ (19) 

where W{ki) is the Fourier transform of W{1). The assumption that the size of the cluster 



is much larger than the typical scales of the fluctuations yields \W{ki 
S{ki) is the Dirac delta function, and therefore K"^^"^ oc k" 



Thus, we find «„ ~ a 



27c6{ki), where 



a 



qj- 



In this section, we have found that, in principle, one can estimate the value of crLN,n 
from analysis of X-ray observations. From the observations one measures o"ln,Sx and asx 
and uses them to infer ULN.n, noting that aq = asx — 0.2. Therefore, one can estimate the 
statistical nature of the intrinsic three dimensional fluctuations from two dimensional X-ray 
observations. 
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Fig. 8. — The average of (Tln.Sx over the 256 synthetic clusters as a function of a, for different 
values of the (3 model power-law index, (3. Symbols correspond to different values of a^ as 
in Figure El Each color shows a different value of (3 (black, red, and blue correspond to 
/? = 1.0,2/3, and 0.5, respectively). Solid, dashed, and dotted lines are fits using equation 
(fT6|) . corresponding to /? = 1.0, 2/3, and 0.5, respectively. The top, middle, and bottom sets 
of three different lines indicate cxLN.n = 0.5,0.3, and 0.1, respectively, as indicated in the 
figure. 

3.3. Potential Systematics 

Using mock observations of isothermal (3 models we found a relation between the intrinsic 
inhomogeneity of the three dimensional cluster gas and the fluctuations in the X-ray surface 
brightness. We turn our attention to the effects of departures from this idealized model. 



3.3.1. f3 Model Power-law Index 

In the above description, we have fiducially assumed the /3 model power-law index 
(3 = 2/3. We investigate two other cases, /3 = 0.5, and /? = 1.0, in Figure [HI where we 
show (Tln.Sx as a function of aq for different cases of [3 (colors). The corresponding fits 
using equation (IT6l) are also shown. Although ctln, Sx tends to increase with increasing /3, the 
change is relatively small (< 10%). 



3.3.2. Temperature Structure 



In the above discussion, we assumed isothermality for the ICM. However, the X-ray 
surface brightness also depends on the underlying cluster temperature structure, including 
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Fig. 9. — The distribution of 6sx for five individual clusters including the effects of temper- 
ature structure. Synthetic clusters (solid histogram) and best-fit lognormal model (dashed 
lines) are both shown for each cluster. Each color corresponds to a different individual syn- 
thetic cluster. Although we display only one example of the power-law index, aq = —3.0, 
similar results are also obtained in other cases. 



a non-isothermal average temperature profile and local inhomogeneity. We investigate these 
effects for the X-ray surface brightness distribution. 



We assume a polytropic profile for the temperature radial distribution expressed as 



T{r) = To 



n(r 



no 



7-1 



(20) 



with polytropic index 7 = 1.2 and Tq = 6 keV, which is the typical set of values in simulated 
clusters (Paper I). The ensemble average of the power spectrum of 6t is assumed to have 
a power-law form ((Pt(^)) ~ Pq{k) oc fc""?-^). Because ar ~ ag,T for the same reasons as 
described in § l3.1l for density fiuctuations, we fiducially adopt the power-law index a^^ = —3 
based on the results of cosmological hydrodynamic simulations (for details see 



We create the lognormal distribution 6t for temperature fiuctuations in the same manner 
as for the density fiuctuations described in § 13.11 The temperature of an arbitrary point is 
assigned according to 

T(r)=5r(r)T(r). (21) 

We adopt ctln.t = 0.3, because it is the typical value for simulated clusters (Paper I). In 
addition, we assume that 6n and 6t are distributed independently, following Paper I. The 
X-ray surface brightness is given by 



5x(R) oc [n{T)]^A[T{r)]dl, 



(22) 
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Fig. 10. — The effect of the PSF on ctln, Sx as a function of radius, -R/rc, for the case of 
Uq = —3 and /3 = 2/3. Sohd curves show crLN,Sx(-R) without convolution of the PSF. Dashed, 
dash-dotted, and dotted curves correspond to 6'hpd/^c = 0.1,0.2 and 0.5, respectively. Two 
values of o"LN,n are plotted, 0.1 and 0.5, as indicated in the figure. 



where A(T) is th e X-ray cooling funct ion. We calculate A(T) in the energy range 0.5-10.0 keV 
using SPEX 2.0 (JKaastra et al.lll996l ) on the assumption of coUisional ionization equilibrium 
and a constant metallicity of 30% solar abundances. 

Examples of the distribution of 6sx in individual clusters are shown in Figure M (solid 
histogram) along with the best fit lognormal distributions (dashed lines). Each color corre- 
sponds to a different individual synthetic cluster. Although only one value for the power-law 
index, Oq = —3, is shown, similar results are obtained for other values. The radial dependence 
of o"ln,Sx including the effects of temperature structure is shown in Figure H] (dot-dashed). 

There are only small differences between the isothermal and non- isothermal cases. The 
X-ray surface brightness depends on the density squared but roughly as a/T for bremsstrahlung 
emission. Therefore, the temperature structure effects on 6sx are much less important than 
those of the density structure. Hereafter, we neglect the effects of temperature structure and 
focus only on the effects of density inhomogeneity. 



3.3.3. Finite Spatial Resolution 

Actual observations by X-ray satellites have finite spatial resolution, characterized by 
the point spread function (PSF). We assume that the PSF is a circularly symmetric Gaussian 
with standard deviation a. The PSF can then be parameterized by a single parameter called 
the half power diameter (6'hpd) in which 50% of the X-rays are enclosed (6'hpd/o" = 2-\/21og2). 
We investigate three cases, ^hpd/6'c = 0.1, 0.2 and 0.5. Figure [TO] shows the effect of the PSF 
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on ctln.Sx as a function of radius. In each case, the average over 256 synthetic clusters is 
shown. Resuhs for no PSF correction (6'hpd = 0, sohd) and 6'hpd/^c = 0.1 (dashed), 0.2 
(dot-dashed), and 0.5 (dotted) are shown. As 6'hpd/6'c increases, o"ln, Sx near the center of 
the cluster decreases. This can be understood as follows. In each radial shell, fluctuations 
smaller than roughly the radius of the shell predominately contribute to the fluctuations, 
namely crLN,Sx(-R)- The PSF effectively smooths out the smaller scale fluctuations (roughly 
up to the size of the PSF), reducing o"ln, Sx, while preserving the large scale fluctuations. 
Since the inner shells only contain small scale fluctuations, they are more strongly affected 
by the PSF. The case of 6'hpd/6'c = 0.5 best illustrates these effects. The reduction of 
(Tln.Sx from the PSF is seen at all radii. However, it is only a slight reduction at large radii, 
increasing as the radius decreases, with a very large effect near the cluster center. 

In summary, when (5„ in three dimensions follows the lognormal distribution, 6^^ in 
two dimensions also approximately follows the lognormal distribution. The mean value of 
cln.Sx for an individual cluster is strongly dependent on both crLN,n and aq. Because aq is 
approximately equal to asx, in principle, one can infer ciLN.n from uln, Sx although there is 
still some dispersion even if aq is known. In addition, the effect of the temperature structure 
is minimal. 



4. Application to Abell 3667 

Simulations suggest that the lognormal model (Eq. [3]) is a reasonable approxima- 
tion of the small scale structure in galaxy clusters. We compare this model with Chan- 
dra X-ray observations of the nearb y galaxy cluster Abell 3667 at a redshift z = 0.056 



( IStruble. M. F.. fc Rood. J. H Ill999l ). A3667 is a well observed nearby bright galaxy cluster 
that does not exhibi t a cool core observed by Chandra . W ith its complex structure, in- 
cludi ng a cold front (IVikhlinin. Markevitch. &: Murray 1 1200 ll ) and possible merger scenario 



Knopp. Henry, fc Briel Ill996l ). A3667 will serve as a difficult test case for the lognormal 



(e.g.. 

model of density fluctuations. 



4.1. Data Reduction 

Chandra observations of the galaxy cluster A3667 are summarized in Table [TJ Listed are 
the observation identification numbers, exposure times, and pointing centers of each of the 
eight archival Chandra observations of A3667 used in this analysis. The data are reduced with 
CIAO version 4.0 and calibration data base version 3.4.2. The data are processed starting 
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Fig. 11. — Chandra image of tlie galaxy cluster Abell 3667 (left) and the corresponding 
^Sx image (right). The counts image has been divided by the exposure map to yield X-ray 
surface brightness (cnt s~^ cm~^ arcmin"^), including scaling for the pixel size. Point sources 
in the field have been masked. 



with the level 1 events data, removing cosmic ray afterglows, correcting for charge transfer 
inefficiency and optical blocking filter contamination, and other standard corrections, in 
addition to generating a customized bad pixel file. The data are filtered for ASCA grades 0, 
2, 3, 4, 6 and status=0 events and the good time interval data provided with the observations 
are applied. Periods of high background count rate are excised using an iterative procedure 
involving creating light curves in background regions with 500 s bins, and excising time 
intervals that are in excess of 4 a from the median background count rate. This sigma 
clipping procedure is iterated until all remaining data lie within 4 a of the median. The 
final events list is limited to energies 0.7-7.0 keV to exclude the low and high energy data 
that are more strongly affected by calibration uncertainties. Finally, the images are binned 
by a factor of eight, resulting in a pixel size of 3.94". This pixel size matches the resolution 
of the synthetic clusters considered in ^ In particular, the ratio of pixel size to the cluster 
core radius of the Chandra image is similar to the synthetic c luster grid spacing cornpared 
to the synthetic cluster core radius, namely, for 6^ ~ 180" (JReiprich fc Bohringer I l2002l : 
Knopp. Henry, fc Brief Ill996l ). 9pi^/9c ~ c?grid/'^c ~ 0.02. Exposure maps are constructed 



for each observation at an energy of 1 keV. The binned images and exposure maps for each 
observation are then combined to make the single image and exposure map used for the 
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analysis. 

A wavelet based source detector is used to find and generate a list of potential point 
sources. The list is examined by eye, removing bogus or suspect detections, and then used 
as the basis for our point source mask. Figure [11] (left) shows the Chandra merged image 
of A3667, the counts image divided by the exposure map, where the point source mask 
has been applied. Also shown is the ^ sx image (right), discussed below. A cold front 
( IVikhlinin. Markevitch. fc Murray II2OOII ) is clearly visible in the south-eastern region of the 



(5sx image. 



4.2. Analysis and Results 

In order to determine the center of A3667, a (3 model is fit to the data with fixed core 
radius (180") and (3 (2/3), using software original ly developed for the combined analysis of 



X-ray and Sunyaev-Zel'dovich effect observations (JReese et al.ll200Cll . l2002l : iBonamente et al. 



20061 ). Because A3667 is nearby and appears very large, Chandra observations do not en- 
compass the entire cluster but provide a wealth of information on the complexities inherent 
in galaxy cluster gas. By using a P model fit to the diffuse emission of the cluster gas we ob- 
tain a better measurement of its center than simply using the brightest pixel or other simple 
estimates, which fail to take into account the complex structure manifest in this cluster. A 
circular region of radius ~ 8' centered on A3667 is used in the analysis, corresponding to two 
and a half times the cluster's core radius, the largest usable region from the arrangement of 
the combined Chandra observations. 

The average X-ray surface brightness is required to compute ^sx = Sx/Sx- If one 
computes the average surface brightness, Sx, in annular shells, then one will tend to under 
(over) estimate Sx toward the inner (outer) radius of each annulus. Therefore, this will 
lead to an over (under) estimate of 6^^ toward the inner (outer) radius of each annulus. To 
alleviate this systematic, we adopt the azimuthally averaged X-ray surface brightness as the 
model for Sx, and use cubic spline interpolation between radial bins. The X-ray surface 
brightness radial profile for A3667 is shown in Figure [121 along with the interpolated model 
(line). 

The probability distribution of 6sx, p(<^Sx), is computed from the histogram of pixels 
calculated from the 5sx image and shown in Figure [T3l The lognormal distribution (Eq. [3]) 
is fit to the p{Ssx) of A3667, where the only free parameter is the standard deviation of the 
logarithm of (5sx, ctln.Sx- The best fit value for the lognormal model is (Tln, Sx = 0.30. In 
addition, a Gaussian distribution is also fit to the data, with its usual two parameters, the 
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Fig. 12. — Chandra radial profile of the galaxy cluster Abell 3667 (points) with the in- 
terpolated model (solid line). This model is used as the average X-ray surface brightness 
distribution in the calculation of 6sx- 
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Fig. 13. — Probability distribution of 6sx from Chandra observations of the galaxy cluster 
Abell 3667 (blue histogram) along with the best fit lognormal distribution (red line) with 
c"ln,Sx = 0.30. The lognormal distribution seems to be a reasonable description of the ICM 
inhomogeneity in A3667. Also shown are the best-fit Gaussian model (dashed green) and a 
Poisson model (dot-dashed magenta) using the average counts per pixel within the fitting 
region. 



mean and standard deviation. Figure [13] shows the PDF of (5sx for the Chandra observa- 
tions of the galaxy cluster A3667 (solid blue histogram). The best fit lognormal (solid red) 
and Gaussian (dashed green) models are also shown. A Poisson distribution (dot-dashed 
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Fig. 14. — Power spectrum of 5sx (thick solid) from Chandra observations of the galaxy 
cluster Abell 3667, normalized to one at the largest scale. Also plotted are three power- 
law power spectra with spectral indices of -2 (dashed), -3 (dot-dashed), and -4 (dotted) for 
comparison. 



magenta) is also shown for comparison, using the average counts per pixel in the fitting 
region as the parameter for the Poisson distribution. Clearly, what is seen is not the result 
of Poisson statistics. The lognormal model seems to be a reasonable match to the observed 
PDF. 

However, without information on the power spectrum of the ^sx fluctuations, it is difficult 
to interpret the value of uln.Sx ( §3.2.2p and relate it to the fluctuations in the density 
distribution (Eqs. [T6l[T7] : Fig. [6]). Therefore, we take the Fourier transform of the 6sx image 
and compute the average power spectrum in wavenumber annuli. The power spectrum of 
6sx fluctuations is shown in Figure [H] (thick solid) along with three power-law spectra with 
spectral indices of -2 (dashed), -3 (dot-dashed), and -4 (dotted) for comparison. The power 
spectrum of ^sx has been normalized to one at the largest scales. A simple power-law model 
fit to the power spectrum yields a spectral index of asx = —2.7 using the entire spectrum, 
and a spectral index of asx = —3.0 if excluding the larger wavenumbers (> 2 arcmin"^), 
roughly where the power spectrum changes shape. 



4.3. Implications 

Both the standard deviation of the logarithm of X-ray surface brightness fluctuations. 



c"LN,Sx = 0.30, and the power spectrum power-law index asx 



-3, fall into the range 
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Fig. 15. — An example of a 6sx map from a cosmological hydrodjTiamic simulated cluster 
( "Centaurus" ) both before (left) and after (right) removal of a quadrant with a large clump. 
Circles show the projected virial radius (-R2oo)- Although within the projected virial radius, 
R200, these structures often reside outside of the three-dimensional virial radius, r2oo- 

expected from hydrodynamical galaxy clusters and therefore used in the synthesized cluster 
analysis ( §2.21) . By combining these pieces of information, we can relate the information 
obtained from the X-ray surface brightness distribution to that of the underlying density 
distribution, using the results of the synthesized cluster analysis. Using the synthetic cluster 
result that the spectral indices of the X-ray surface brightness fluctuations and that of the 
Gaussian field are simply related by asx ~ a? + 0.2, and the relation between o"LN,n, ctln, Sx, 
and ttq (Eqs. [121 [H]; Fig. E]), the Chandra results of ctln.Sx = 0.30 and asx = —2.7 imply 
that the fluctuations in the underlying density distribution have aLN.n = 0.43. A value 
of asx = —3.0 implies o"LN,n = 0.36. The difficult test case of the A3667 X-ray surface 
brightness seems to follow the lognormal distribution of density fluctuations, thus enabling 
an estimate of the statistical properties of the underlying ICM density fluctuations. 



5. Application to the Cosmological Hydrodynamic Simulated Clusters 



Results from cosmological hydrodynamic simulations motivated the lognormal model 
for ICM inhomogeneity. In ^ we found that synthetic clusters with lognormal fluctuations 
show a linear relation between uln.Sx and cxLN.n- We now return to clusters extracted from 
cosmological hydrodynamic simulations to further explore these results. 

For each cluster extracted from the simulations, we create X-ray surface brightness 
maps towards three orthogonal directions, and compute 5sx(R) = Sx(R)/Sx{R) in a similar 
manner as described for the synthetic clusters in § I3.1.2[ The regions we consider are within 
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Fig. 16. — The distribution of 5sx for each of the six clusters from a cosmofogical hydrody- 
namic simulation (points and solid histogram). Each color indicates the projection along a 
different, orthogonal line of sight. For each line of sight, we show the number of quadrants 
used for the analysis. For example, "3/4" indicates that one quadrant is excluded and three 
remain. The best fit lognormal model for each projection is also shown (dotted lines). 



the projected virial radius R2oo- The projected virial radius, -R20O) is the radius within which 
the mean interior density is 200 times that of the critical density. 

Although the lognormal distribution is a good fit to the density (and temperature) of 
simulated galaxy clusters in three-dimensions, the projection to X-ray surface brightness 
suffers from the additional complexity of projection effects. If large clumps are present, the 
distribution of X-ray surface brightness fluctuations, ^sx, is not well approximated by the 
lognormal distribution. The large clumps artificially distort the average profile of the cluster 
and therefore bias the value of 5sx, which depends on the average profile. We also note that 
although these clumps fall within the projected virial radius, -R2005 they usually fall outside 
of the three dimensional virial radius, r2oo- We therefore exclude quadrants that contain 
large clumps, using 6sx > 10 as the exclusion criterion. Then, we recompute S'x(-R) and 
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Fig. 17. — The density fluctuation standard deviation predicted by our model, 
(TLN.n (model) = crLN,Sx/<5(a(j) versus that from the simulations, (jLN,n(sim). Symbols show 
different simulated clusters (see figure legend) and colors indicate different orthogonal fines 
of sigfit. Also plotted is tfie simple linear relation crLN,n(model) = (jLN,n(sini) for comparison. 



^Sx- Tfie complex structure of simulated clusters is illustrated in tfie 5sx images sfiown in 
Figure [151 wfiere examples of a simulated cluster botfi before and after removal of a quadrant 
are displayed. Tfie circles sfiow tfie projected virial radius, -R2oo- 

In Figure [16] tfie probability distributions of 6sx for tfie simulated clusters (fiistograms) 
along witfi tfie best-flt lognormal model (dotted lines) are displayed. Eacfi color indicates tfie 
projection along a different, ortfiogonal line of sigfit. Overall, tfie probability distributions of 
(5sx are reasonably well approximated by tfie lognormal function, consistent witfi tfie results 
from tfie syntfietic clusters (§ l3.2.ip . 

We now come full circle to compare our results from tfie syntfietic clusters directly 
to tfie simulations. In order to do tfiis, we look at tfie relationsfiip between (TLN,n(sim) 
measured in tfie simulated clusters and o"LN,n (model) predicted from tfie syntfietic cluster 
results, equations ([T6|) and ([T7|l . wfiere we adopt a^ = asx — 0.2 (see §3.2.2p . Tfie value of 
asx for eacfi simulated cluster is obtained by fitting a power-law model, P{K) oc 7^"^", to 
tfie power spectra of (5sx- Because tfie resolution of tfie simulations is mucfi poorer tfian tfiat 
of tfie syntfietic clusters, we must recompute tfie coefficients Ci and C2 in equation [17] from 
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a set of lower resolution synthetic clusters. Assuming Vc ~ 100 h^^ kpc for the simulated 
clusters, we choose the resolution ~ O.lrfgrid/rc, noting that this value corresponds to the 
maximum resolution of the simulations. Performing the same procedure described in §31 we 
obtain ci = 3.99 x lO'^ and Cs = 3.36 x lO'^. 

We compare ciLN.n (model) and crLN,n(sim) in Figure [T71 Each color corresponds to a 
different line of sight. Although there is large scatter, these results indicate that it is possible 
to estimate cxLN.n within a factor of two only using the information obtained from the X-ray 
surface brightness distribution. 



6. Summary 

We have developed a method of extracting statistical information on the ICM inhomo- 
geneity from X-ray observations of galaxy clusters. With a lognormal model for the fluc- 
tuations motivated by cosmological hydrodynamic simulations, we have created synthetic 
clusters, and have found that their X-ray surface brightness fluctuations retain the lognor- 
mal nature. In addition, the result that ctln.Sx and cxLN.n are linearly related implies that 
one can, in principle, estimate the statistical properties of the three dimensional density 
inhomogeneity (aLN.n) from X-ray observations of galaxy clusters (ctln.Sx and asx). 

We have compared the predictions of our model to Chandra X-ray observations of the 
galaxy cluster A3667. For the first time in a real galaxy cluster we were able to detect the 
lognormal signature of X-ray surface brightness fluctuations, which was originally motivated 
by simulations. Based on the synthetic cluster results, this enabled an estimate of the 
statistical properties of the inhomogeneity of the ICM of A3667. In the context of lognormally 
distributed inhomogeneity, we obtain crLN,n ~ 0.4 for the gas density fluctuations of A3667. 
It is encouraging that the value of the fluctuation amplitude for Abell 3667 is in reasonable 
agreement with typical values from the simulated clusters. 

Finally we check the validity and limitation of our method using several clusters from 
cosmological hydrodynamic simulations. Unlike the fairly idealized synthetic clusters, simu- 
lated clusters exhibit complex structure more akin to real galaxy clusters. As a result, the 
empirical relation between the two- and three-dimensional fluctuation properties calibrated 
with synthetic clusters when applied to simulated clusters shows large scatter. Nevertheless 
we are able to reproduce the true value of the fluctuation amplitude of simulated clusters 
within a factor of two from their two-dimensional X-ray surface brightness alone. 

Our current methodology combined with existing observational data is useful in de- 
scribing and inferring the statistical properties of the three dimensional inhomogeneity in 
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galaxy clusters. The fluctuations in the ICM have several implications in properly interpret- 
ing galaxy cluster data. In particular, our current model may be useful in interpreting data 
from current and future galaxy cluster surveys using the Sunyaev-Zel'dovich effect, which 
have the potential to provide tight constraints on cosmology. 
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A. Relation of Density and Surface Brightness Distributions Under the 

Thick-slice Approximation 

Modeling galaxy clusters with a spherical isothermal (3 model (Eq. [TT]) . the surface 
brightness at an arbitrary projected angular radius, 9, is given by 



s^{e) oc 



[n{v)fdl 



M2n^Q 1 + 



4^2 \ -3/3 .oo 



5nnU) 



1 + 



rc,esi^) 



-3/3 



dL 



(Al) 



where rc,eff(6') = ^/rf+d^9^, and we assume the ctlx.x in equation [3] is independent of r. 
Therefore, the second moment of n (M2 = (n^)/(n)^ = exp (— ctln „)) is also independent of 
r. In the above, we use (n) = n{r). 



The ensemble average of {Sx{0)) can be expressed as 

J —00 



^ 1.. r(3/?-l/2) / d\B 



•2q2\ -3/3+1/2 



(A2) 
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Combining equations (lAll) and flA2l) . (5sx,ens reduces to 



^Sx,ens(6') = Kf3 (5nn(r) 



rcM^) 



-3/3 



rf 



TcM^) 



K/3 



TT 



-1/2. 



r(3/5) 



r(3/5-i/2)- 



(A3) 



Now, fixing 6*, let us consider the three-dimensional field 5„„(r)iya(/) and its projected 
two-dimensional field i^sxensle, defined as 



(^Sx.ensle 






(A4) 
(A5) 



where we use a dimensionless length normalized by rc,eflf(^) distinguished by prime (/' 
l/rc,es{0), k'l = kiTc^esiQ))- Then, we can consider the variance of the 5sx,ens|e-field, 



a 



'Jsx.GnsIs 



— i-^ I dK'Ps^ l(K'), 



where Psy^^^^\g{K') is the (two-dimensional) power spectrum of (5sx,( 
5nn field can also be written as 

1 



a 



Onn 



(2vr) 



rfk'P„„(k'). 



(A6) 
The variance of the 

(A7) 



With this, the relation between Ps^^^^(K') and Pnn(k-') is 



P5x,..|.(K') = - / dk[PU^')\W,{k[)\\ 



(A8) 



The Fourier conjugate Wp^kl) is given by 



w^i^'i) = '^M ^ 



-3/3+1/2 



27? 



K^ 



3/3+1/2 



(H) 



(A9) 



.,/ r(3/3)- 

where i^_ 3/3+1/2 (/i;|) is modified Bessel function of the second kind . 

In the case that the largest s cale fiuctuation is smaller than the physical scale (the 
thick-slice approximation, following iFischera fc Dopita I (120041 )). the Fourier conjugate of the 
window function becomes the Dirac-delta function, |VF^(A;9P ~ 9{l^)^{^'i)- The normalization 
factor g{f3) is given by 

-r(3/5)2r(6/? - 1/2) 



»""-V„ *'l"'><*''l' = ^^r(3/3-i/2pr(6^)- 



(AlO) 



-30- 
Let us define the effective width 



Aeff(^) = 27rr,,eff(^)/^(/5) = y^^!^y^/^^^f<^MO). (Ali: 



Fischera fc Dopita I (l2004j ) explore the column density distribution assuming a plane parallel 
geometry with width A. In the thick slice case, Acg corresponds to A although they consider 
the column density not the surface brightness. We assume statistical isotropy and a power 
law spectrum with upper and lower limit {k'^i,-^ = fcmax''^c,eflf(6') and fc^^j^ = A;min'^c,eflf(6')), 

^ 1=0 Otherwise. ^ ' 

Finally, using equation flA6l) . flA7l) . and flASP under the thick-slice approximation, we 
obtain 

lK„ + 3)(l-C""+2)/Aefr(^)^^' , ^ . , ^ 

OLnn T -^ and ann T -^ 



2 / 2 

<5sx,cns|e/ <5,i 



''max 



2K„ + 2)(l-C"""+=^) V h 



"max 



2(C - 1) V ^: 
1-1/C Meff(^)V' 



U-r)rj ^1 



2 log C V ^max 

where C = /cmax/^min and /max = 27rA;~i^j^. 

Then, although (y\ ^\ is the variance of ^sx.ensle-field, we regard it as the variance 
of the ensemble average of ^sx.enslR-') at Q. The conversion to the standard deviation of 
logarithm is expressed as 



^LN,Sx = Jlog(l + aL.„,J (A14) 



<5sx,cns|e 



In Figure HJ we adopt C, = -^box/(2(igrid) = /s//Ny, where /s and /nj are the sampling 
frequency and the Nyquist frequency, respectively, and /max = -^box- We also adopt the fitted 
value of ann in equation (JAISP . 



B. Distribution of the Density Squared 

If one assumes density inhomogeneity fluctuations, 5„ = n/{n), follow the lognormal 
distribution, PLNlf^n; crLN,n), the fluctuations of the density squared, 6nn = n^ / {'n?) can be 
written as 

5nn = 5n^-ri^ = 5n^eicg,{-al^n)- (Bl) 
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where ( ) indicates ensemble average. 

For simplicity, we replace 5n and 5nn by x and y, respectively, 

x = 5n]y = 5nn- (B2) 

Therefore, the relation between x and y is 

x = y/yexp{(Tl^j2). (B3) 

Because x follows PhN{x; aLN.n), the distribution of y is obtained by 

(ix 

P{y) = PLt^{x; CrLN,n)3- = Pln(2/; 2crLN,n)- (B4) 

dy 
Therefore, Snn follows the lognormal distribution with lognormal standard deviation of 

2crLN,n- 
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Table 1. A3667 Chandra Observations 



t(,xp RA DEC 

obsID (ks) (h m s) (dm s) 



513 


45 


20 12 50.30 


-56 50 56.99 


889 


51 


20 11 50.00 


-56 45 34.00 


5751 


131 


20 13 07.25 


-56 53 24.00 


5752 


61 


20 13 07.25 


-56 53 24.00 


5753 


105 


20 13 07.25 


-56 53 24.00 


6292 


47 


20 13 07.25 


-56 53 24.00 


6295 


50 


20 13 07.25 


-56 53 24.00 


6296 


50 


20 13 07.25 


-56 53 24.00 



